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The staggered 6- vertex model describes the competition between surface roughening 
and reconstruction in (100) facets of CsCl type crystals. Its phase diagram does not 
have the expected generic structure, due to the presence of a fully-packed loop-gas 
line. We prove that the reconstruction and roughening transitions, cannot cross nor 
merge with this loop-gas line if these degrees of freedom interact weakly. However, 
our numerical finite size scaling analysis shows that the two critical lines merge 
along the loop-gas line, with strong-coupling scaling properties. The central charge 
is much larger than 1.5 and roughening takes place at a surface roughness much 
larger than the conventional universal value. It seems that additional fluctuations 
become critical simultaneously. 
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1. INTRODUCTION 



In a recent paper, Mazzeo, Carlon , and van Beijeren [1] discuss the competition 
between surface roughening and reconstruction in c(2x2) reconstructed (100) facets 
of CsCl type crystals. Their numerical finite size scaling (FSS) results for the 
staggered 6-vertex model are quite surprising. The phase diagram lacks a so-called 
reconstructed rough (RR) phase, although such a phase is a generic feature in 
surfaces where step excitations do not destroy the reconstruction order. The phase 
diagram should have the same structure as for missing row (MR) reconstructed 
simple-cubic (SC) (110) facets [2-4]. The roughening and reconstruction lines should 
be able to cross. Instead they only seem to approach each other exponentially 
close. In this paper, we explain why this happens. The absence of a RR phase in 
the staggered 6-vertex model is accidental, the result of a special symmetry of the 
interactions in this particular model, the presence of a fully-packed loop-gas line. 

In section 2, we review the rich history of the staggered 6-vertex model. In sec- 
tion 3, we describe the topological properties of step and wall excitations in c(2x2) 
reconstructed CsCl(lOO). We set up the cell-spin model description for this type of 
competition between surface roughening and reconstruction. Topological consider- 
ations determine whether the roughening and reconstruction lines can cross or only 
merge (whether a RR phase is possible or not). For example, in MR reconstructed 
SC(110) facets they can cross, but in MR reconstructed FCC(llO) facets they can 
only merge [2-4]. We show that in c(2x2) reconstructed CsCl(lOO) they are allowed 
to cross. The competition in this surface is in the same universality class as in MR 
reconstructed SC(110) facets. However, the RR phase in the staggered 6-vertex 
model is at best narrow. We estimate the energies of the two topologically distinct 
types of steps, and find that they cost almost the same energy in the region of the 
phase diagram where the surface roughens. 

Carlon et al. [5,6] point out the existence of a special line in the phase diagram. 
It runs in-between the Ising and roughening lines. Along this line the partition 
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function is equivalent to the 4-state Potts model on a square lattice with negative 
Boltzmann weights. They expect this line to be similar to a so-called disorder line 
and that this will explain their numerical results (the non-crossing of the Ising and 
roughening lines) . In section 4 we show that this line is equivalent to a fully-packed 
loop-gas on a square lattice. 

In section 5, we prove rigorously that the reconstruction line cannot cross 
the loop-gas line. Furthermore, we show that the roughening line cannot cross 
the loop-gas line either if the roughening and reconstruction degrees of freedom 
couple weakly. The weak-coupling hypothesis assumes that the reconstruction and 
roughening degrees of freedom interact weakly, such that their scaling properties 
are a superposition. Earlier studies of models for MR reconstructed SC(llO) and 
FCC(llO) facets strongly support the weak-coupling hypothesis [2-4]. It should 
hold for CsCl(lOO) as well since the cell-spin model of section 3 is the same. 

This seems to resolve the issue. The phase diagram found by Mazzeo et al. 
[1] is the only one allowed within weak-coupling theory, but it is an accident. The 
special symmetries of the loop-gas line cause this particular model to follow a special 
cut through the generic phase diagram. The roughening and Ising lines approach 
each other only pathologically close, because entropy cannot be lowered far enough 
to reach the crossing point into the RR phase. In general, the interactions in CsCl 
type surfaces will be more generic and allow the RR phase. However, this is not 
the end of the story. 

In section 6 we present our numerical FSS results. Mazeo et al. [1] performed 
their study before they discovered the loop-gas line. Knowledge of the exact location 
of the line where the roughening and reconstruction lines must merge or cross (if they 
do so) enhances the accuracy of the analysis considerably. We find that the scaling 
behaviour along the loop-gas line does not obey the weak-coupling hypothesis. 

The question of weak versus strong coupling type competition between recon- 
struction and roughening degrees of freedom is an important unresolved issue in the 
theory of 2D critical phenomena. It appears not only in surface physics, but also 
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e.g. in coupled Josephson junction arrays in a magnetic field (the fully frustrated 
XY model) [7-13]. The phase diagrams of these problems share as basic feature a 
conventional order-disorder transition line (such as an Ising or 3-state Potts transi- 
tion), approaching a critical (rough) phase. A critical line described by conformal 
field theory (CFT) with central charge c < 1 competes with a critical phase with 
central charge c = 1. The fundamental question is whether new interesting CFT's 
can result from this competition. The c-theorem [14,15] implies that such CFT's 
must have a central charge larger than one. The ones we know are rather simple 
direct products of c = 1 and c < 1 theories, some with extra symmetries, such 
as super-symmetry, between the two types of degrees of freedom. This supports 
the weak-coupling hypothesis. However, do more interesting types of c > 1 CFT's 
really not exist? Can the coupling between two types of c < 1 degrees of freedom 
ever lead to more intricate c > 1 scaling behaviour? This question is the driving 
force behind a large number of numerical studies, in particular within the context 
of the fully frustrated XY model and the competition between surface roughening 
and reconstruction [1-13]. But it proved to be extremely difficult to answer. 

Three types of answers are possible and each has appeared in the literature. 
The first one is that the c < 1 line cannot reach the c = 1 phase. They approach each 
other only pathologically close. This suggests the existence of a no-go theorem of 
some sort. Our exact results in section 5 amount to such a no-go theorem, but only 
for the staggered 6-vertex model, and only within the weak-coupling hypothesis. 

The second possibility is the weak-coupling scenario. The reconstruction and 
roughening lines cross or merge, but the two types of critical fluctuations interact 
weakly and behave like a direct product. The central charge is equal to the sum 
(c = 1.5 in our case, since c = 0.5 at Ising critical points, and c = 1 inside the rough 
phase) . This type of behaviour is almost indistinguishable from the effective scaling 
when the two lines approach each other only pathologically close (see section 5). 
Typically, the numerical data can be interpreted both ways [7-13]. Fortunately, our 
exact results in section 5 distinguish between the two in the 6-vertex model. 
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The third possibility is that the lines cross or merge with scaling behaviour 
different from a simple superposition. The central charge is not equal to the sum. 
Convincing evidence for strong-coupling would revolutionize CFT at c > 1. Numer- 
ical evidence for strong-coupling has been presented in models related to the fully 
frustrated XY model, but remains ambiguous [7-13]. We find strong numerical 
evidence (in section 6) that the roughening and reconstruction lines in the stag- 
gered 6-vertex model merge along the loop-gas line, with strong-coupling scaling 
properties. 

In section 7, we summarize our results, discuss related recent results for fully- 
packed loop-gases on different lattices, and give a possible explanation for the strong 
coupling behaviour. 



5 



2. THE STAGGERED 6- VERTEX MODEL 



The (100) facets of CsCl type crystals have a body centered type stacking with 
two kinds of atoms, type A and B. The appropriate solid-on-solid description is a 
staggered Body-Centered Solid-on-Solid (BCSOS) model, equivalent to a staggered 
six-vertex model. Stacks of A atoms occupy the A-sublattice where the column 
heights are odd numbers, Ha = ±1, ±3, ±5, .... The B type atoms occupy the B- 
sublattice where the column heights are even, h B = 0, ±2, ±4, .... Nearest neighbour 
stacks differ in height by only one unit, dh = ±1. The surface energy is given by 

W = { E A( h (x,y) ~ /l( x +l,y+l)) 2 + E B (h( x+ljy ) - h( X)y+1 )) 2 }. (1) 

x,y 

The summation runs over only the A- type sublattice sites (all even values of x + y). 

This model has a rich history. Fig.l shows its phase diagram. It is exactly 
soluble by the Bethe Ansatz along the line Ea = E B , where it reduces to the 
so-called F-model [16-18]. We denote the flat phase as F(n + §), because the 
average surface height is an half-integer. A more detailed notation is in terms of 
the heights at 4 sublattices, (fiA,+, h>A,-i hs,+, hs,-) = (n + 1, n + 1; n, n). The + 
and — indices represent the two (checkerboard type) sublattices within each A and 
B sublattice. Along Ea = Eb, elementary step excitations induce a height change 
±1. They roughen the surface at za = zb = |, with za = exp(£ , ^//cBT) and 
z B = exp(E B /k B T) [16-18]. 

Ea and E B are not equal in CsCl type surfaces, since the A and B type atoms 
interact differently. Knops [19] realized that this changes the disordering of the 
F{n + |) phase. He used the equivalence of eq.(l) to the Ashkin- Teller model. 
Rephrased from the more recent perspective of preroughening (PR) transitions and 
disordered flat (DOF) phases, his results are as follows: the F(n + ^) contains two 
types of steps, Sa and Sb- In surfaces with a height 2m+ \, up-steps are of type Sb 
and down-steps are of type Sa (they reverse roles in surfaces with a height 2m — \). 
—Ea is the energy of Sa steps (per unit length), and — Eb of Sb steps. The free 
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energy of Sa steps vanishes before that of Sb steps on the Ea > E B side of the 
phase diagram. This does not cause roughening yet, since building-up a slope in 
the surface requires Sb steps as well. Only the spontaneous symmetry breaking 
between pairs of surface heights, 2m + | and 2m — |, is being lifted. The result is a 
disordered flat (DOF) phase with lots of Sa type steps in the surface, but still flat 
at large length scales. In Fig.l, we denote this as the DOF(2m) phase, because the 
average surface height is an even integer, 2m. 

This is an example of a preroughening (PR) transition, but in a different uni- 
versality class than the conventional one [2,4]. In both cases the average surface 
height changes spontaneously by half a unit. At conventional PR transitions the 
number of degenerate equivalent surface heights does not change. They all shift by 
one half unit, and the distance between them remains the same. In this example 
however, the distance between degenerate equivalent surface heights increases from 
one to two. The elementary step height is equal to dh = ±1 in the F(n+|) phase, 
but equal to dh = ±2 in the DOF(2m) phase. Therefore, this is a simple Ising line 
instead of a line with continuously varying critical exponents. 

This doubling in the basic step height not only changes the nature of the transi- 
tion into the DOF phase, it also delays the roughening transition considerably, even 
for small energy differences Ea — Eb [19]. The surface roughness is characterized 
by the amplitude of the height-height correlations as 

((h r+ro -h r ^))^-±-log(\r\). (2) 

Rough phases become unstable towards discreteness of the surface height at Kq = 
7r/8 for step size dh = ±2 compared to Kq = tt/2 for step size dh = ±1, not until the 
surface is four times as rough. In this particular model, these roughening lines lie at 
the side of the phase diagram where both step energies are negative (both Ea and 
Eb are positive). Point F in Fig.l is located at za = zb = y/{l + ±y/2) = 1.30656. 
The local structure of the phase diagram around point F is known as a "critical 
fan" [20]. Notice that unreconstructed surfaces described by eq.(l) never roughen 
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(except along Ea = E B ). They follow specific paths through Fig.l, approximately 
lines at constant Ea/E b with both Ea and E B negative. Such lines do not enter 
the critical fan. The absence of a roughening transition is not a generic feature, 
however. Experimental unreconstructed CsCl(lOO) type surfaces will include step- 
step interactions and other aspects that are able to move point F towards and 
beyond za = z B = 1. 

In upper left corner of Fig.l the surface reconstructs. All B-type columns are 
at the same height, 2m, and the A- type columns alternate between 2m ± 1; e.g. 
(hA+,h.A-,h B+ ,h B -) = (2m + 1,2m — 1,2m, 2m). We call this the RA(2m,9) 
reconstructed phase. The average surface height is an even integer. The Ising type 
order parameter 9 = 0, it denotes which of the two A-type sublattices is on top. 
The average surface height is the same as in the DOF(2m) phase. The difference is 
the appearance of anti-ferromagnetic type ordering of the A-sublattice heights. The 
competition between the deconstruction of this type of order and surface roughening 
is the topic of this paper, in particular the existence of point S in Fig.l, and the 
scaling properties of the critical line beyond this point. 

Fig.l has mirror symmetry with respect to Ea = E B . The A and B type 
particles switch roles. In the lower right corner of Fig.l the surface height is an odd 
integer. The surface reconstructs into the R B {2m + 1, 9) phase and the DOF phase 
is of type DOF(2m+ 1). 
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3. STEPS AND WALLS IN THE C(2x2) PHASE 



Reconstructed surfaces can disorder in several ways: they can lose their recon- 
struction first and then roughen; the can roughen first and only later deconstruct; 
or roughening can induce a simultaneous deconstruction transition. The latter hap- 
pens when the topology of the surface implies that steps destroy the reconstruction 
order parameter. Fig. 2 shows the three topologically distinct line excitations in the 
i?A(2m, 9) phase: a (+2, 7r) step, a (0, n) wall, and a (—2,0) step. Wall excitations 
do not change the surface height. They cause a switch in which of the two A- type 
sublattices is on top. Steps of type (±2, 0) change the surface height by ±2 but do 
not change the Ising order parameter 6>; the same A-type sublattice stays on top. 
Steps of type (±2, n) change the surface height by ±2 and switch which of the two 
A-type sublattices is on top. 

It is an illusion to think that (±2, 7r) steps destroy the reconstruction order. 
They preserve a different definition of the reconstruction order. It can be expressed 
in terms of which sublattice is on top, 9 = 0, tt, or in terms of parity spins, Sa = 
exp^m/i^] = ±1. The Sa spins are ordered anti-ferromagnetically in the i?yi(2m, 9) 
phase. These two definitions of the Ising order are equivalent in flat surfaces but 
inequivalent in rough surfaces. The (±2, tt) steps destroy the sublattice order but 
preserve the parity order. The (±2,0) steps destroy the parity order but preserve 
the sublattice order. Two types of reconstructed rough (RR) phases are possible: 
the surface is rough, but such that the Sa order parameter remains non-zero; or 
the surface is rough, but such that the 9 order parameter persists. In diffraction 
experiments, the roughening transition can be easily mistaken for a simultaneous 
deconstruction transition. The "reconstruction diffraction peak" couples only to 
one of the two order parameters (typically the parity one) and that order might not 
be the type that persists in the RR phase. 

Along paths where either type of step costs less energy than a wall, the surface 
roughens first, into the appropriate RR phase, followed by a Ising type deconstruc- 
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tion transition inside the rough phase. Along paths where walls cost less than steps, 
the surface deconstructs into the DOF(2m) phase, and only later roughens. 

Fig. 3 illustrates these different typical behaviours. It is a schematic phase 
diagram for the following model [4]. Each site of a square lattice contains an Ising 
spin, cr r = ±1, and an height variable , h r = 0, ±2, ±4, ... . They interact as 

H = - Yl K ° a r(?r> + (K T a r a r ' + K h ) [2 - (h r - /v) 2 ] (3) 

<r,r'> 

with r and r' nearest neigbour sites. Only steps of height dh = ±2 are allowed. 
cr r = cos(7r# r ) = ±1 represents the sublattice type reconstruction order parameter 
and h r the local surface height. Walls cost E(0,n) = 2{K a + K T ). Steps cost 
respectively E(2, 0) = 2(K T +K h ) and E(2, n) = 2(K a +K h ). This is a renormalized 
effective model, on a larger length scale than the staggered 6-vertex model. It must 
be in the same universality class in the local neighbourhood about the reconstructed 
phase, assuming we identified correctly all critical fluctuations of the 6-vertex model 
in this part of the phase diagram. Fig.l should follow a specific 2D cut through 
Fig.3. 

The mirror symmetry in Fig.3 with respect to A = E(2, 0) — E(2, n). reflects 
the equivalence between the two definitions of the reconstruction order. Consider 
the following construction of a typical configuration. Define a second type of Ising 
spin S r = exp(^iirh r ) = ±1, and draw cr and S type Bloch walls along the bonds of 
the lattice. The a type Bloch walls represent wall excitations in the surface and the 
S type Bloch walls represent (2, 0) type steps. Place one arrow along each S type 
Bloch wall to denote the direction in which the height changes across the steps (up 
from left-to-right while looking in the direction along the arrow). Sections along 
the lattice where S and a type Bloch walls coincide represent (2, tt) type steps. The 
following change of variables generates the mirror symmetry in Fig.3. Define new 
Ising spins r = a S to represent the (2, ir) type steps and then eliminate the a spins, 
(a, S) — > (r, S). This leads us back to eq.(3), but with the r spins replacing the a 
spins, and with K a <-> K T . The two RR phases switched place. 
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This transformation is reminiscent of super symmetry, between the fermion 
(Ising) and boson (height) degrees of freedom. It is weaker than super symmetry, 
only a Z2 type invariance [4]. At the same time it is more general, an exact sym- 
metry of the lattice model not restricted to T c . The entire A = space is invariant, 
instead of only the critical line A-B. We expect that critical fluctuations generate 
full super symmetry at large length scales. The roughening induced simultanous de- 
construction transition along A-B will then be described by a super symmetric CFT, 
probably one with central charge c = 1.5 where the roughening and reconstruction 
are weakly coupled. 

Experimental systems and microscopic models follow specific cuts through 
Fig. 3. For example, the anti-ferromagnetic RSOS model describes checkerboard 
type reconstructed SC(IOO) facets [2,4]. Indeed, its phase diagram represents a 
generic slice out of Fig. 3 with a DOF phase and a RR phase; and the Ising and 
roughening degrees couple weakly with central charge c = 1.5 [2,4]. One of the 
exactly soluble generalized RSOS models [21] moves along the Ising surface in Fig. 3 
as well and confirms weak-coupling behaviour [22]. A third example is the chiral 
4-state clock-step model [4] which describes MR reconstructed FCC(llO) facets. 
Topology requires the two types of steps in those surfaces to have identical ener- 
gies. The non-chiral limit of the 4 state clock-step model coincides with the A = 
plane of eq.(3). Numerical evidence supports the expectation that along A-B the 
Ising and roughening degrees of freedom couple weakly with c = 1.5. 

The phase diagram of the staggered 6-vertex model should to be a generic 
cut through Fig. 3 similar to the RSOS model. There is no intrinsic topological 
requirement for A to be zero. However, Mazzeo et al. do not find a RR phase. 

One possible explanation is that A is small or vanishes "accidentally" in the 
6-vertex model. The wall and steps in Fig. 2 run diagonally across the surface. In 
that direction, the two types of steps cost the same amount of energy per unit 
length, £(2,0) = E(2,ic) = \V2 (2E A - E B ). A wall cost E(0,ir) = E A . 
This suggests that A is equal to zero. However, A is quite large for walls and steps 
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running in the horizontal and vertical direction: (0, 7r) walls cost E(0, n) = Ea 
per unit length, (2, 0) steps cost E{2, 0) = Ea — E B , and (2, it) steps exist only as 
composite objects (a (2, 0) step followed by a (0, 7r)-wall, i.e., E(2, n) = IE a — Eb)- 

Walls prefer to run in the diagonal direction, but steps switch direction. Deep 
inside the R(2m, 9) phase to the left of the line Ea/Eb — —0.7 the steps prefer to 
run in the diagonal direction. A is small, but this is the part of the phase diagram 
where the walls cost much less energy. The surface deconstructs into the DOF(2m) 
phase before it roughens. Along z B = the model reduces to the Ising model, and 
therefore the deconstruction transition takes place at za = 1 + V2- In zeroth order 
approximation, the deconstruction line is located at za = 1 + \/2 for all Eb since 
the wall energy does not involve E B - 

Deep inside the R{2m, 9) phase to the right of the line Ea/E b — —0.7 steps pre- 
fer the vertical and horizontal directions. The (2, 0) steps are most favourable, and 
A is large. However, roughening cannot take place until Eb and Ea are of the same 
order of magnitude. Roughening takes place at approximately exp[£'(2, 0)//cbTr] 
~ 1 + \[2 (the Ising formula gives reasonable estimates for transition temperatures 
in general). We can construct two different estimates for the roughening line in 
Fig.l, by assuming the (2,0) steps run vertically or diagonally. These estimates 
are quite close to each other. This means that near the roughening transition the 
E(2, 0) steps run almost equally likely in the diagonal as in the horizontal or vertical 
directions. A must be small near the roughening transition. E(2, it) steps come into 
play, and the RR phase is narrow at best. 



12 



4. FULLY-PACKED LOOP-GAS ON A SQUARE LATTICE 



Carlon et a/.[5,6] realized recently that along the lines za + zb = 1 and za = 
zb ± 1 (the dashed lines in Fig.l) the staggered 6-vertex model maps onto the 
4-state Potts model and that the za = zb + 1 line seems to lie in-between the 
Ising and KT roughening lines. This mapping has been known actually for a long 
time but not from this perspective. For details, we refer to the original source [23]. 
Carlon and van Beijeren [6] expect that the Potts line will turn out to be a type of 
disorder line, and thus will explain that the Ising and roughening lines cannot meet, 
in accordance with their numerical results [1] . The properties of this line are much 
more intriguing. The essential observation is that along the Potts line the 6-vertex 
model reduces to a fully-packed (FP) loop-gas. 

In the 6-vertex representation of the BCSOS model, an arrow points along 
each bond of the lattice, to denote the height difference between nearest neighbour 
columns Iia — h B = ±1. Fig. 4a shows the 6 allowed vertex states. In the loop- 
gas model every bond contains a loop segment. The loops follow the bonds of the 
lattice, are closed, and do not intersect. It is a fully-packed loop-gas. Fig. 4b. shows 
the 2 possible vertex states, A and B. The partition function is of the form 

z = J24 a 4 b i Nl = z% v J2 xNb2Nl i ( 4 ) 

Q Q 
with the summation over all FP loop-graphs Q, x = zb/za, Na the number of 
vertices of type A, Nb the number of vertices of type B, Ny = Na + Nb the 
total number of vertices in the lattice, and N L the number of loops. The fugacity 
factors of 2 can be counted by placing arrows on the loops; clockwise and anti- 
clockwise arrows. Loop configurations with such arrows resemble configurations in 
the 6-vertex model, but they are not one-to-one related. The vertex states with 
anti-parallel arrows, labeled 5 and 6 in Fig. 4a, can be interpreted as both A or 
B type loop states. How to deal with this is one of the essential steps of the 
mapping of the Potts model onto the 6-vertex model: each 6-vertex configuration 
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represents the sum over all possible loop-gas interpretations [23]. Still, the models 
are only identical along special lines. In the 6-vertex model, vertex states 5 and 
6 are assigned a Boltzmann factor uj^ = ujq = 1, (both next nearest neighbour 
heights are equal), while in the loop-gas their weights are the sum over all loop 
interpretations: u$ = ujq = za + zb • The staggered 6-vertex model reduces to the 
FP loop-gas, eq.(4), only when these are equal, only along the line za + zb = 1. 

All lines of type ±za + ±zb = 1 are FP loop gases as well, due to the fact 
that vertex states 1 and 2 (and also 3 and 4) in Fig. 4a appear always in pairs (Fig.l 
has mirror symmetry with respect to za and also zb)- The lines za = zb + 1 and 
zb = za + 1 represent loop-gases with negative Boltzmann weights 

z = J2(-i) Nb 4 a 4 b z Nl = 4 V J2(~ 1 ) NbxNb2Nl - ( 5 ) 

g g 

These two lines are analytic continuations of each other with < x < 1 along 
za = zb + 1 and 1 < x < oo along zb = za + 1; the minus signs in eq.(5) can be 
counted equally well by Nb as Na since Na + Nb = Ny is a constant. FP loop- 
gases have been a focus of attention recently [24-26]. In particular, the FP loop gas 
on a honeycomb lattice resembles eq.(4). We will discuss possible connections with 
this recent work in section 7. 

In the loop-gas, the arrows are merely a gimmick to count the loop fugacity. 
They are placed at random on each loop. Therefore it seems reasonable that any 
order associated with the up-down nature of the steps must be absent along the 
loop-gas line, the reconstruction order as well as the surface flatness. This is too 
naive. The surface is able to maintain flatness. Only the reconstruction order is 
absent. Notice that the loop gas lines in Fig.l move through the DOF phases. The 
following visualization is quite useful. Interpret the B sublattice as "patches of red- 
land" and the A sublattice as "patches of blue-land" . The loops are the coastlines. 
The two vertex states in Fig. 4b represent the presence of either a bridge connecting 
the two patches of blue-land or the two patches of red-land. The fugacity factor 2 
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for each loop can be interpreted as a random height difference of dh = ±1 between 
blue- and red-lands at each coastline; while looking in the direction along the arrow 
the land on the left is lower by one unit than the land on the right. (The arrows 
attribute a helicity to each coast line.) At za = there exist only red-land bridges. 
All red patches are connected and they are all at the same height; the surface stays 
flat. All blue patches are disconnected and randomly distributed at heights h±l. 
Consequently, the reconstruction order is absent. Everywhere along the loop-gas 
line inside the DOF(2m) phase, there exists a large continent of red-land spanning 
the entire lattice, and keeping the surface flat. Moreover, all blue-lands are finite 
in size (lakes inside the red-land continent). The red- and blue-lands switch roles 
inside the DOF(2m + l) phase. The only other possibility is that all red and all blue 
land masses are finite in size. There the surface is rough and unreconstructed. The 
reconstruction order is always absent along the loop-gas line. In the next section 
we prove this rigorously. 
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5. INTERFACE FREE ENERGIES 



Consider the 6-vertex model partition function in a semi-infinite strip geometry, 
and the following boundary conditions (BC's): h(x + N,y) = ±h(x,y) + a, with 
a an even integer. The lattice forms a cylinder, infinitely long in the y-direction 
and N lattice sites in circumfence in the x-direction. The free energy per unit strip 
width for each of these BC's, /(±,a), can be calculated as the largest eigenvalue 
of the transfer matrix. The free energy differences, r/(±, a) = /(±, a) — /(+, 0), are 
related to specific step, wall, and other defect free energies. 

For example, the h(x + N,y) = h(x, y) + 2 BC forces a step into the _R/i(2m) 
reconstruction. To see this is, vizualize the RA(2m) ground state as a criss-cross 
structure of horizontal and vertical intersecting straight lines. Each carries an arrow, 
pointing alternately up and down (to the left and right). Reversing all arrows 
interchanges the two degenerate i?^(2m) states. Reversing the arrows along only 
one vertical line creates a step. This is a (2, 0) type step, not a (2, tv) type step, 
because the arrow reversal along this specific line does not affect the 6>-type order 
parameter on either side of the step. The same A-type sublattice stays on top on 
either side. By reversing the arrow we create a net height difference across the 
surface of a = ±2. The h(x + N,y) = h(x, y) + 2 BC matches this structure for 
even values of N, and therefore forces a (2, 0) step into the surface along the entire 
cylinder. //(+, 2) is equal to the (2, 0) type step free energy everywhere inside the 
i?^(2m) phase. 

One method to force a wall excitation into the reconstructed phase is to apply 
periodic boundary conditions (PBC), h(x + N,y) = h(x,y), at odd strip widths 
N. To see this is, visualize the i2^(2m) ground state as an array of elementary 
loops with alternate helicity. In the red/blue-land interpretation of the loop-gas 
the i?A(2m) ground state is the structure, in which all blue-lands are disconnected 
elementary lakes, and the height changes at the coast lines follow a strict up-down 
pattern. The coast line arrows have alternate helicity. We run the transfer matrix 
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in the diagonal direction, where the reconstructed phase fits only onto the lattice 
if the strip width N is a multiple of 2. The h(x + N,y) = h(x, y) BC frustrates 
the helicity order at odd strip widths. Therefore r/(+, 0) o is equal to the wall free 
energy. 

The BC h(x + N, y) = h(x, y) + 2 forces a (2, n) type step into the reconstructed 
phase at odd strip widths. Other types of BC's have similar effects: twist boundary 
conditions (TBC) at even values of N create a (0, n) wall for h(x + N, y) = —h(x, y), 
and a (2, n) type step for h(x + N,y) = —h(x, y) + 2. 

Loop-gas lines in solid-on-solid models signal special properties. Free energies 
with certain boundary conditions become "accidentally" equal, implying that spe- 
cific excitations have identical free energies. For example, the RSOS model contains 
a (non fully-packed) loop-gas line, which coincides with the exact location of the 
roughening line. Its presence proves the existence of the preroughening transition 
in that model [2]. 

The partition function of the loop-gas does not change when we modify the rules 
for placing the arrows on the loops. Consider the TBC, h(x + N,y) = —h(x 7 y). 
The seam is the vertical line across the entire cylinder where this boundary condi- 
tion is being implemented. (Its location is gauge invariant; moving the seam and 
deforming it does not alter the partition function.) The columns on one side of the 
seam interpret the columns on the other side as being at height —h. In the arrow 
representation, this means that the direction of the arrow on each loop reverses 
each time that loop crosses the seam. There are two types of loops, homotopic 
and non-homotopic loops. Non-homotopic loops wrap around the cylinder in such 
a way that they cannot be contracted topologically into a point (like beads on a 
necklace). The requirement that the arrow on each loop reverses each time it crosses 
the seam is incompatible with non-homotopic loops. The configurations with PBC, 
h(x+N, y) = h(x, y), are almost identical to those with TBC, h(x+N, y) = —h(x, y). 
Each homotopic PBC configuration is matched by a TBC configuration. Moreover, 
their Boltzmann weights are identical, since along the loop-gas line the reversal 
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of arrows at the seam does not affect the Boltzmann factor. However, all non- 
homotopic PBC configurations are absent for TBC. Therefore, the free energies 
/(+, 0) and /(— , 0) are identical except for the contribution to /(+, 0) of configura- 
tions with non-homotopic loops. Such configurations are suppressed in the RaC^tti) 
and DOF(2m) phase because all loops are finite (all blue-lands are finite sized lakes). 
This means that rj(— , 0) = /(+, 0, ) — /(— , 0) vanishes in the thermodynamic limit 
exponentially fast with N. The latter is incompatible with reconstructional order, 
since r}{— , 0) represents the free energy of a wall excitation in the Ra^tti) phase, 
and therefore cannot be equal to zero. The loop-gas cannot lie inside the i?yi(2m) 
phase. 

To make this argument rigorous, we add the following aspect to the boundary 
conditions. Draw a second seam and associate a phase factor </> = n/2 ((f) = — 7r/2) 
each time a loop crosses this seam with an arrow pointing from left to right (right 
to left) . These phase factors do not affect the Boltzmann factor of homotopic loops, 
(the phases add-up to zero), but they freeze-out non-homotopic loops (the phases 
add-up to \n mod(7r)). Define free energies ?y(±, a, 4>) = /(±, a, </>) — /(+, 0, 0), with 
a = 0, ±2, ±4, and = 0, ir/2. From the above discussion it follows that along the 
loop-gas line 

t/(- 0,0) = 7/(+,0,^7T) (6) 

for all even strip widths N. The twist boundary condition on the left creates a 
frustration of the Ising order, while the periodic type boundary condition on the 
right is compatible with Ising order. This proves rigorously the absence of long- 
range Ising order along the loop-gas line. The loop-gas line cannot enter the Ra(2iti) 
reconstructed phase nor either of the two RR phases. It cannot cross the Ising line 
in Fig.l. 

This is a rather weak result. In particular, it does not imply the absence 
of a RR phase because it does not exclude the possibility that the loop-gas line 
enters the rough phase. The implications of eq.(6) become much stronger within 
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the constraints of the weak-coupling hypothesis. Next, we summarise the scaling 
properties along every path the loop-gas can follow through Fig. 3 and confront 
each with the loop-gas symmetry eq.(6). This summary is important also for the 
numerical analysis in section 6. 

(i) Suppose the loop-gas line enters the deconstructed rough phase. The central 
charge is equal to c = 1, and all the above surface free energies decay at large N as 
the inverse of the strip width with amplitudes that are linked to each other as: 



The inverse of Kg, defined in eq.(2), represents the surface roughness. The rough 
phase is described by the Gaussian fluctuations at large length scales, and these 
relations are exact and easy to derive in the Gaussian model. Kq must be smaller 
than Kq < 7r/8 since the KT roughening transition takes place at Kg = n/8. 
This is a factor 4 smaller than the conventional value Kg = tt/2, because the step 
excitations create a height difference dh = ±2 instead of dh = ±1. The loop-gas 
identity eq.(6), applied to eq.(7), yields the value Kg = 7r/2, inconsistent with 
Kg < 7r/8. The loop-gas line cannot enter the deconstructed rough phase. (This 
argument proves also that point P in Fig.l must coincide with the KT transition 
into the F(n + ^) phase.) 

(ii) Suppose the loop-gas line moves along the Ising surface inside the rough 
phase, in particular the phase boundary with the RR phase dominated by (2, 0) 
type steps. At Ising critical points the central charge is equal to c = 0.5 and the 
Ising Bloch wall free energy scales as a power law with universal amplitude 



N Vs (+,a, 4>) 




Nrj s (-,a, 0) 



TV 



(7) 



4 



Nrji ~ 7r/4. 



(8) 



The central charges and universal amplitudes add-up, as c = 1.0 + 0.5 = 1.5, when 
the Ising and roughening degrees of freedom interact weakly. The roughness degrees 
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of freedom behave like in eq.(7), with Kg < 7r/8, and the Ising degrees of freedom 
like in eq.(8). This leads to the following FSS amplitudes 



N V (+, 2,0) ~Nr) s {+,2,0)~2K G 



iV77(- 0,0) ~iV[7/ a (-,0,0) + 7/i] 



1 



iVr/(-,2,0) ~iV[7/ a (-0,0) + 7/i] ~ -7T 
A^(+,0,0) o -ATr/i ~ ^tt 

iVr/(+, 2, 0) o ~AT[r/ s (+, 2, 0) o + m] ~ 2K G + ^tt. 




(9) 



The loop-gas identity eq.(6), applied to eq.(9), yields the value Kg = 7r/4, still too 
large compared to Kg < 7r/8. The loop-gas line cannot move along the Ising plane 
inside the rough phase. Along the opposite Ising plane, the phase boundary with 
the RR phase dominated by (2, n) steps, the scaling relations are similar, with the 
two types of steps reversing roles. 

(iii) Suppose the loop-gas line follows the line segment A-B in Fig. 3. The 
roughening and Ising degrees of freedom still couple weakly, The central charge 
remains equal to c = 1.5, and the wall free energy scales still as Nrji ~ ^7r, but the 
surface roughenss is constant, Kg = 7r/8 [4]. The FSS amplitudes are similar to 
those in eq.(9). The major difference is that the amplitudes for the two types of 
steps are identical by symmetry, 



The two interface free energies in eq.(6) behave the same as in (ii). This yields 
again the value Kq = 7r/4, inconsistent with Kq = tt/8. The loop-gas line cannot 
move along line segemnt A-B. 

(iv) Suppose the entire loop-gas line lies inside the DOF phase. The FSS 
central charge estimates decay to c — > 0. Both step free energies are non-zero. The 
Ising wall type free energy is equal to zero. However, the convergence becomes 



N V (+2, 0) ~ N V (+, 2, 0) o ~ 2K G . 



(10) 
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extremely slow when the Ising and roughening lines approach each other closely. 
The asymptotic forms will not be reached, and the apparent scaling will be almost 
indistinguishable from (ii) or (iii). For A^O, the effective scaling behaviour will be 
almost identical to that in eq.(9), but still be distinguishable; the surface roughness 
cannot exceed the value K g = tt/8. For A = 0, the effective scaling behaviour will 
be identical to that for (iii) . This is a fundamental dilemma in studies of this type 
of phenomena. The only distinction between scenarios (iii) and (iv) is a judgement 
call on whether the two lines merge or not in the numerical analysis. In our case, 
eq.(6) resolves the issue. It excludes (ii) but allows (iv). 

In summary, the loop-gas suymmetries exclude all the above scenarios except 
the last one. The only possibilities are: either the roughening and Ising lines ap- 
proach each other pathologically close; or the roughening and reconstruction degrees 
of freedom interact strongly and thus circumvent the above arguments. In any case, 
the Ising line can never cross the loop-gas line, since that aspect does not require 
the weak-coupling hypothesis. 
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6. Numerical Results 



To distinguish between the scenarios outlined in section 5, we calculate the 
exact free energies /(±, a, </>) for semi-infinite lattices of width 2 < N < 10, for 
several boundary conditions as defined in section 4. We run the transfer matrix in 
the diagonal direction, where the state vector is 2 2N dimensional. Such strip widths 
are in the usual range for transfer matrix calculations. They might seem small for 
readers more familiar with Monte Carlo simulations, but realize that our values of 
/(±, a, 4>) are accurate to better than 12 decimal places. There is no statistical noise, 
unlike MC simulations. This allows a very detailed finite size scaling (FSS) analysis 
that incorporates the leading corrections to scaling. It pays to trade system size 
for the ability to determine the leading corrections to scaling, because at criticality 
FSS corrections decay only as power laws. We known the exact location of the line 
where the reconstruction and roughening transitions must merge or cross (if they 
do), the FP loop-gas line za = z B + 1. This makes our numerical analysis more 
accurate than earlier studies of the same type of phenomena. 

Fig. 5 shows the FSS estimates for the central charge c(N) along the loop-gas 
line, za = zb + 1- These values follow from the free energy with periodic boundary 
conditions /(+, 0, 0) and the conformal field theory scaling relation 

f(+,0,0) N ~f o +^c (11) 

as 

c(N) = ^ (iV2 4 ~ 1)2 [/(+, 0, OV-i - /(+, 0, 0) N+l ]. (12) 

Eq.(ll) is valid at criticality. The c(N) approximants of eq.(12), must converge 
to zero away from criticality (exponentially fast). Indeed they do so inside the 
DOF(2m) phase at small za- At criticality c(N) must converge to the value char- 
acteristic for the universality class of the phase transition. In Fig. 5, we present the 
raw c(N) data together with our best FSS estimates for c (the dashed line) and 
rather conservative error bars (the shaded area). At large za, the FSS corrections 
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to scaling become large and the FSS analysis of the type c = c(N) +A/N X becomes 
less reliable. The essential point is that c is certainly larger than c = 1.5. This 
contradicts all weak-coupling scenarios (see section 5). Maybe c varies continuously 
along the loop-gas line, but it is more likely that c is a constant between 2 < c < 3 
and that the variation of c in Fig. 5 reflects crossover scaling behaviour between that 
value and c = 1.5 at the point where the roughening and Ising lines merge. 

Fig. 6 shows the FSS scaling behaviour of Nrj(+, 0, 0) o and Ni](— , 0, 0) along 
the loop-gas line. The first one forces an Ising wall into the reconstruction, the 
second one an Ising wall and a twist in the surface. Both free energies vanish in the 
DOF(2m) phase as expected. At large za they scale with amplitudes that converge 
towards the values Nrj(+, 0, 0) o ~ n/A and Nr}(— , 0, 0) ~ 7r/2 (the dashed lines) 
consistent with the weak-coupling scenarios (ii)-(iv). 

Fig. 7 shows the FSS behaviour of Nrj(+, 2, 0) and Nrj(+, 2, 0) o along the loop- 
gas line. Nrj(+, 2, 0) induces a (2, 0) type step in the surface and Nrj(+, 2, 0) o a (2, n) 
type step. Both diverge in the DOF(2m) phase as expected (the step free energies 
are finite). At large za, both decay as power laws. These data are inconsistent 
with scenarios (i) and (iii), in which Nr](+ : 2, 0) and Nrj(+,2,0) o should become 
equal at large N. The results are inconsistent with scenario (ii) as well, since 
Nr)(+, 2, 0) and Nrj(+,2,0) o should differ by tt/4 when the loop-gas line moves 
along the Ising plane inside the rough phase (see eq.(9)). To demonstrate this, 
we plot in Fig.7 N[rj(+, 2, 0) o - rj(+, 0, 0) o ] instead of Nrj(+, 2, 0) o itself. The two 
sets of curves should fall on top of each other and converge towards a continuously 
varying roughness parameter 2Kq- Instead, they differ by a factor close to 2. Most 
importantly, the data is inconsistent with scenario (iv), in which the roughening and 
Ising lines only approach each other asymptotically close. We should find effective 
scaling behaviour of the form Nr](+, 2, 0) ~ N[r](+, 2, 0) o - r/(+, 0, 0) o ] ~ 2K G with 
an effective roughness Kq > ir/8. Not only do the two amplitudes differ by a factor 
2, both sink well below 2Kq = tx/A. The surface roughness becomes too large by a 
factor of about two. 
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Fig. 8 shows an example, of the FSS behaviour of Nrj(+, 2, 0) along the cut 
za + zb = 11 through the loop-gas line. There exists no RR phase on either side 
of the loop-gas line. In the reconstructed phase, Nrj(+, 2, 0) must diverge and in 
the rough and RR phases converge to 2Kq . At the KT roughening transition, 
2Kq must be equal to 2Kq = 7r/4. The conventional method to determine the 
roughening temperature is to extrapolate the points where Nr)(+, 2,0) = n/4 as 
function of N. In Fig. 8, these points lie at the reconstructed side of the loop-gas 
line. They converge towards the loop-gas line at such a rate that power law fits 
actually overshoot the loop-gas line. This might lead to the conclusion that the 
roughening line does not cross nor merge with the loop-gas line. On the other 
hand, along the loop-gas line itself, the amplitude converges very well to a value 
much smaller than 7r/4. This behaviour is similar to what happens at a conventional 
KT type roughening transition if one tries to estimate T c by extrapolating the points 
where the FSS amplitude is larger than the true critical value, Nrj(+, 2, 0) = n/4 + a 
with a > 0. The FSS behaviour in Fig. 9 strongly supports the absence of a RR 
phase. Nr](+, 2, 0) diverges everywhere on the reconstruction side of the loop-gas 
line. Roughening seem to take place exactly at the loop-gas line but surprisingly at 
a surface roughness well above the universal KT value. 

Fig. 9 shows two types of estimates for the critical point where the rough- 
ening and Ising lines merge: the Nrj(— , 0, 0) crossing points from Fig. 6, and the 
Nr}(+, 2, 0) = 7r/4 points from Fig. 7. The existence of crossing points in Nr](— 7 0, 0) 
and Nr)(+, 0, 0) o (see Fig. 6) is quite significant. In the more global context of Fig.l, 
these crossing points converge to the Ising type reconstruction lines. (This is the 
conventional method to locate such critical lines.) Their presence along the loop- 
gas line implies that all FSS estimates for the Ising transition cross the loop-gas 
line. The existence of points along the loop-gas line where A^?y(2,0, 0) = 7r/4 (see 
Fig. 7) is equally significant. These points converge to the KT roughening lines in 
Fig.l. (This is the conventional method to determine such roughening lines.) Their 
presence along the loop-gas line implies that all FSS estimates for the roughening 
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transition cross the loop-gas line as well (see also the discussion about Fig. 8). The 
Ising and roughening lines cross the loop-gas line from opposite directions. All finite 
N estimates for the Ising and roughening lines cross each other. 

Both curves in Fig. 9 must converge to E c — > oo if the roughening and Ising 
lines only approach each other asymptotically close. This is very unlikely, although 
both curves are convex (the corrections to scaling behave as Ea(N) = E c + A/N x 
with an exponent x < 1). A conservative estimate puts the critical point somewhere 
between 1.4 < E c < 1.6. 

These numerical results contradict all weak-coupling scenarios, in particular the 
one where the roughening and Ising lines only approach each other pathologically 
close, the only one allowed by the loop-gas symmetries of section 5. The most 
damaging evidence is that Nr](+, 2, 0) becomes too small by a factor 2. It seems 
too far fetched that this amplitude can rebound all the way back to tt/4 at very 
large N. Moreover, the central charge is significantly larger than the weak-coupling 
value c = 1.5. Finally, the FSS estimates for the Ising and roughening line cross 
each other for all finite N, and lead to an estimate E c ~ 1.5 ± 0.1 for the point 
along the loop-gas line where they merge. 
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7. CONCLUSIONS. 



In this paper, we study the phase diagram of the staggered 6-vertex model from 
the perspective of the competition between surface roughening and reconstruction, 
and also in the context of unresolved issues about the scaling properties of QFT 
with central charge c > 1. In section 3, we review the weak-coupling hypothesis 
as encountered in previous studies for this type of interplay. In particular we show 
that a reconstructed rough (RR) phase must be present in typical c(2x2) CsCl(lOO) 
type surfaces. The staggered 6-vertex model includes a fully-packed (FP) loop-gas 
line (section 4). Its special symmetries explain the absence of the RR phase in this 
particular model. In section 5, we prove that the Ising type reconstruction line 
cannot cross the loop-gas line. Moreover, we demonstrate that within the context 
of the weak-coupling hypothesis the roughening line cannot cross the loop-gas line 
either. This would explain the results by Mazzeo et al. and put them in agreement 
with the generic phase diagram Fig. 3. However, our numerical results contradict all 
weak-coupling scenarios. The Ising and roughening lines merge along the loop-gas 
line. The central charge is large, at least equal to c = 2 (see Fig. 5), and the surface 
roughness increases to a value about twice as large as the universal KT value (see 
Fig.7). 

Eqs.(4) and (5) represent FP loop-gases on a square lattice. The FP loop-gas 
model on a honeycomb (HC) lattice and a four-coloring problem on the square lattice 
are related to this. Those models have large central charges as well, respectively 
c = 2 and c = 3 [24-26]. Unfortunately they seem more closely related to eq.(4) 
than eq.(5). The partition function of the FP loop-gas on a HC lattice is similar 
to eq.(4), with three types of "bridge-energies" instead of two, z a with a=A,B,C 
(the plaquettes of the HC lattice form three sublattices instead of two). Point 
P in our phase diagram is a critical point with central charge c = 1. It is the 
meeting point of two DOF phases and also the point where the Ea = Eb surface 
roughens (see Fig.l). The corresponding point in the FP loop-gas on the HC lattice, 
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za = zb = zc , is a critical point with central charge c = 2 [24-26]. It is the meeting 
point of three DOF phases, and of three critical lines with central charge c = 1 (the 
phase boundaries between pairs of DOF phases) [27]. The scaling properties along 
za = zb ± 1 in Fig.l appear to be more complex than this. This is not surprising 
since eq.(5) includes negative Boltzmann weights. 

We close with some speculations about the origin of strong-coupling scaling 
along the loop-gas lines za = zb ± 1. It might represent a novel type of conformal 
field theory with c > 1. More likely, it represents a conventional CFT, but one 
in which more degrees of freedom become critical than in eq.(3) and Fig. 3. The 
obvious candidate is the Ising degree of freedom of the R(2m + 1, 9) phase located 
on the opposite side of the Ea = E b line in Fig.l. The two reconstructed phases 
do not seem to meet in Fig.l, but they actually do so via the back-door. They 
meet at point x = 1 in eq.(5), since za = zb + 1 and zb = za + 1 are each other's 
analytic continuations. A 4-state clock model coupled to a SOS model describes 
both types of reconstructions simultaneously. The a = ±1 Ising spins in eq.(3) 
denote which of the two A-type sublattices is on top in the R(2m, 9) phase. Their 
generalizations are 4-state clock variables, 9 = 0, ±^7r, 7r, that point in the direction 
of the polarization of the arrows in vertex states 1-4 in Fig. 4a. They denote which 
of the four sublattices is no top: one of the two A-type sublattices, 9 = 0, it; 
or one of the two B-type sublattices 9 = ±^7r. Such a model has ample room for 
conventional CFT's with central charge c > 2. The loop-gas symmetries will enforce 
a non-generic path through its phase diagram. Point x = 1 in eq.(5), where the two 
reconstructions meet, is almost certainly a critical point with central charge c = 2. 
Fig. 3 applies when x = 1 is an isolated critical point. In that case, the Ising and 
roughening lines cannot meet until x = 1 (our exact results of section 5). Instead, 
our numerical analysis shows that x = 1 is not an isolated point. The next step 
will be an analytic calculation, to determine the scaling properties at point x = 1. 
The critical dimension of the crossover operator in the loop-gas direction must be 
irrelevant or marginal for the loop-gas to remain critical until point £ in Fig.l. This 
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line segment is probably some some sort of Baxter-line coupled to roughening, since 
the universal amplitudes of both step free energies in Fig. 7 vary continuously. 
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FIGURE CAPTIONS 



1. Phase diagram of the staggered 6- vertex model defined in eq.(l). At point S 
the Ising reconstruction and KT roughening lines merge into a single transition 
along the loop-gas line (dashed). 

2. Excitations inside the i?^(2m, 9) reconstructed phase: a (2, 0) step, (0, tv) wall, 
and (2, tt) step running in the diagonal direction. 

3. Generic phase diagram for the competition between surface roughening and 
Ising type reconstruction, see eq.(3), with A the energy difference between 
(2, 0) and (2, n) type steps, and R the difference between the energy of a (0, tv) 
wall and the average step energy. 

4. The six vertex states of the 6- vertex model (a), and the two vertex states of 
the loop gas model (b), with their Boltzmann weights. 

5. FSS estimates for the central charge along the loop gas line for strip widths 
N < 10. The shaded area represents a conservative estimate for the uncertainty 
in the extrapolated values (the dashed line). 

6. Ni](+,0, 0) o (free energy of a (0, tc) wall) and Nrj(— , 0, 0) (free energy of a 
(0, 7r) wall with a twist in the surface) along the loop gas line for strip widths 
N < 10. 

7. Nrj(+, 2, 0) (free energy of a (2, 0) step) and Nrj(+, 2, 0) o - Nrj(+, 2, 0) o (free 
energy of a (2, tv) step minus that of a (0, tv) wall), along the loop gas line for 
strip widths N < 10. 

8. The step free energy Nr](+, 2, 0) along the line za + zb = H for strip widths 
N < 10. 

9. FSS estimates for the location of the critical point exp(— E c ) = za = zb + 1 
along the loop gas line from: (a) the Nrj(— , 0, 0) crossing points in Fig. 6, and 
(b) the Nr](+, 2, 0) = tt/4 points in Fig.7. 
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Davidson and den Nijs, Fig. 5 




Davidson and den Nijs, Fig. 6 




Davidson and den Nijs, Fig. 7 




Davidson and den Nijs, Fig. 8 




Davidson and den Nijs, Fig. 9 



